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ABSTRACT 

The force of gravity acting within the volume occupied by young, compact and 
massive superstar clusters, is here shown to drive in situ all the matter deposited 
by winds and supernovae into several generations of star formation. These events 
are promoted by radiative cooling which drains the thermal energy of the ejected 
gas causing its accumulation to then rapidly exceed the gravitational instability 
criterion. A detailed account of the integrated ionizing radiation and mechanical 
luminosity as a function of time is here shown to lead to a new stationary solu- 
tion. In this, the mass deposition rate M, instead of causing a wind as in the 
adiabatic solution, turns into a positive feedback star-forming mode equal to the 
star formation rate. Some of the implications of this extreme positive feedback 
mode are discussed. 

Subject headings: HII regions: clusters: winds - galaxies 

1. Introduction 

The discovery by HST and large ground-based telescopes, of a large population of su- 
perstar clusters (SSCs) in a large variety of galaxies (see Ho 1997, Martin Hernandez et al. 
2005, and Melo et al. 2005), has driven the community to speculate about the impact that 
these new units of star formation may have on the ISM of their host galaxies. SSCs with 
a mass, Msc, ranging from a few xlO 5 M to up to 6 xlO 7 M (see Walcher et al. 2004; 
Pasquali et al. 2004) all within a small volume of radius 3 to 10 pc, are indeed some of 
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the most energetic entities within galaxies. Their large UV photon output and mechanical 
energy are now believed to be the largest negative feedback agents in starburst galaxies, 
leading not only to a large-scale structuring of the ISM and to limit star formation, but also 
to be the agents capable of establishing as in M82 a supergalactic wind, thereby removing 
processed material from galaxies and causing the contamination of the IGM (Tenorio-Tagle 
et al. 2003). 

In the present series of papers however, we have shown that there is a threshold limit 
that massive and energetic compact clusters may cross to find out that radiative cooling 
inhibits their stationary outflow condition (Silich et al., 2004) and then the matter ejected 
by the stellar sources, unable to escape, is to accumulate within the star cluster volume 
(Tenorio-Tagle et al. 2005). The facts in such cases are that radiative cooling drastically 
diminishes the sound speed and the pressure gradient within the thermalized ejected matter, 
inhibiting the possibility of a wind. Radiative cooling upsets the stationary condition in 
which the deposited matter (M) has to equal, at all times, the amount of matter streaming 
out of the SC volume, M = ILscIv 2 ^ = 4:irR sc pscCsc, where Csc is the amount of energy 
injected by stellar winds and SNe, csc is the local sound speed at the SC surface and 
is the resultant wind terminal speed in the absence of radiative cooling. As soon as this 
happens, the ejecta begins to accumulate, promoting larger densities and an even faster 
cooling within the SSC volume. 

The location of the threshold line in the Csc vs size (Rsc) diagram (see Figure 1), the 
line that defines whether or not a wind is inhibited, depends on several variables. It depends 
on the size of the star-forming region (Rsc), the metallicity of the ejected gas, which has a 
strong impact on the cooling curve (see paper III). It also depends on the assumed Csc/M or 
adiabatic terminal speed (v oo) of the wind. The latter is also bound to the usual assumption 
that the energy deposited by SN is always 10 51 erg but the mass of potential SNe within the 
cluster ranges from, say, 100 M Q to 8 M and so the injection speed (similar to v^) and the 
deposited amount of matter are also functions of time. Another factor that strongly affects 
the location of the threshold line is the thermalization efficiency (e) which simply defines 
the fraction of the mechanical energy that can be evenly spread within the cluster volume. 
Estimates of e by several authors lead to values between 1 (Chevalier & Clegg; 1985) to 0.03 
(see Melioli & Del Pino 2004 and references therein) and depends simply on the proximity 
of the sources undergoing winds and SNe, which through radiation may reduce the amount 
of energy available after thermalization. We have shown that there are three different types 
of solutions: SSCs far away from the threshold line (low mass, low energy clusters) undergo 
a quasi-adiabatic evolution well described by the Chevalier & Clegg (1985) and Canto et al., 
(2000) solution. More energetic clusters are to have strongly radiative winds. Cooling hardly 
affects their velocity (v w ~ Voo) and density distribution (p w oc r~ 2 ), but their temperature 
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instead of falling as r~ 4 / 3 , it falls rapidly to T w ~ 10 4 K close to the SC boundary and 
the more so, the closer they are to the threshold line. The strongly radiative winds around 
such clusters lead, compared to the adiabatic solution, to very much reduced X-ray envelope 
sizes. The third solution is for clusters above the threshold line. These would have their 
winds inhibited and as shown below, this turns them into very efficient positive feedback 
star-forming agents. 

Here we follow the evolution of a massive and compact cluster and show the fate of the 
mass reinserted back, through winds and SNe, into the SC volume. Section 2 and 3 give a 
full description of SSC as positive feedback agents and section 4 summarizes our conclusions. 



Here we follow the evolution of a stellar cluster, assumed to be the result of an instanta- 
neous or coeval burst of star formation. The cluster has an initial total mass in stars (M$c) 
equal to 10 7 M within an Rsc = 5 pc. We assume a standard Salpeter IMF and an upper 
and lower mass limits equal to 100 M Q and 1 M , respectively. The mechanical energy 
of the cluster Csc ~ 10 41 erg s^ 1 places our example cluster in the zone where stationary 
superwinds are inhibited (see Figure 1) and thus the mass returned by the stars (M) begins 
immediately to accumulate within the star cluster volume, favoring an even faster cooling. 
Figure 2 shows how the density of the accumulating gas (p ac = 3Mt/AnR 3 sc ; where t the 
evolution time), grows as a function of time within the SSC volume. 

Following this trend, clearly there is a moment when the density of the accumulating 
gas will exceed the gravitational instability criteria 



where p,j is the Jeans density and k,G,p and 7 are the Boltzmann constant, gravitational 
constant, mean mass per particle and the ratio of specific heats (5/3), respectively and Rsci 
is the SSC radius in pc units. At that moment, when p ac = pj collapse will inevitably 
proceed. 

The initial ample supply of UV photons (N° = 10 54 photons s _1 , see Figure 3a), exceeds 
at first the number of recombinations within the volume occupied by the reinserted gas and 
the resultant HII region, given the large metallicities of the ejecta, is here assumed to rapidly 
approach an equilibrium temperature Thii ~ a few 10 3 K. At these temperatures the sound 
speed (< 10 km s" 1 ) remains well below the escape speed (V esc = (2GM S c / Rs) 1 ^ 2 ~ 130 
km s" 1 ) and the reinserted gas would inevitably continue to accumulate to rapidly (within 
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1.5 x 10 6 yr) reach the value of the Jeans density for a gas at say, T HII = 5000 K, and collapse 
into a new stellar generation. The event gives rise to a new phase of matter accumulation, 
which once more will rapidly approach pj (for T— 5000 K) and undergo collapse within 
a free-fall time, of the order of 10 5 yr, while transforming ~ 2 x 10 5 M into stars. All 
stellar generations resultant from mass accumulation within the SSC volume have here been 
assumed to also acquire a Salpeter IMF with similar upper and lower mass limits as those 
imposed to the main superstar cluster, and their resultant properties (mechanical energy 
and UV photon output) have been added to those produced by the main cluster. 

A few, almost identical, stellar generations are expected from the accumulation process 
(solid rising lines in Figure 2), every time that the accumulated gas density p ac reaches 
pj(5000fT) (dashed line in Figure 2). The situation changes slightly when the number 
of ionizing photons (N°), despite the added contribution of secondary stellar generations, 
becomes insufficient to fully ionize the accumulated matter within the star cluster volume. 
This is due to the evolution of the main cluster, whose UV photon output begins to fall 
as t~ 5 after ~ 3.5 Myr (dotted line in Figure 3a). Figure 2 shows pnn (thin solid line), 
the maximum density within the accumulating volume that can be supported fully ionized 
by the UV radiation produced by the evolving cluster (phii — (3iV p? / AttR^qP) 1 / 2 ; where 
p = lArriH and j3, the recombination coefficient to all levels above the ground level = 2.59 
xlO -13 cm -3 s _1 ). During the accumulation process, once p ac exceeds pmi, the ionized 
volume begins to shrink to end up as a collection of ultra compact HII regions around 
the most massive stars left within the cluster, while the bulk of the ejected material, now 
recombined, continues to cool, approaching rapidly (within a time-scale of less than 1000 
yr) a temperature ~ 100 K. Matter is at all times uniformly replenished within the whole 
SSC volume, and thus the gas density presents an almost uniform value. However, the 
accumulating gas now has two different temperatures {Thii an d 100 K) and as p ac grows 
and the fraction of the ionized volume (Jhii — 3A^°/i 2 /47ri?| c /9p^ c ) shrinks, the size of cold 
condensations (at 100 K) able to become gravitationally unstable and their free-fall time 
also become smaller. 

The drop in the number of ionizing photons (Figure 3a) and the consequent growth 
of the neutral volume (Figure 3b) lead then to a second important condition in which the 
characteristic accumulation time r ac = p gas R 3 sc {\ — fmi)/M becomes equal to the free- 
fall time tff — (3ir/32Gp gas ) 1 / 2 . This condition defines p gas , the density above the Jeans 
instability limit for a neutral condensation at 100 K (pj(100K)): 



and thus once p ac (see Figure 2) becomes equal to p gas , a new stationary solution becomes 
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possible. 

Everything happens very rapidly, compared to the evolution time-scale of the parental 
cluster (~ 40 Myr), and almost at the same time. The ejected matter is thermalized within 
the SSC volume and immediately begins to cool. At the same time that it accumulates 
making cooling even faster. This allows it to rapidly reach the required p gas value, above the 
Jeans instability limit, that warrants its collapse in a similar time-scale, while the collapsing 
material is replenished by the newly ejected matter. When this happens, the mass deposition 
rate from the cluster becomes equal to the rate of star formation. Gravitational collapse and 
star formation within the star cluster volume and with the matter injected by all sources, 
drive in this way a new stationary condition through a new era of quasi-continuous star 
formation in which M is equal now to the star formation rate (SFR). 



3. The quasi-continuous star formation era 

All new generations of stars are due to deposit their ejecta from winds and SNe during 
m 40 Myr, the expected duration of the type II SN phase in coeval star clusters (Leitherer 
& Heckman 1995), and at a rate that depends only on the mass of the stellar generation. 
Re-processing leads to an escalating M, as this is directly proportional to the total mass in 
stars t 

M stars (t) = M sc + [ SFR(t)dt (3) 
Jt c 

created within the star cluster volume; where tc marks the start of the quasi-continuous 
star-forming phase and M$c accounts for the initial mass in stars and for all generations 
prior to the quasi-continuous star formation phase (see Figure 2). This promotes a larger star 
formation rate, and a faster fulfillment of the collapse conditions (p ga s > Pj, and T ac — tff). 
Under these conditions the SFR is: 

SFR(t) = ^^ = 5M stars (t)- (4) 
This equation has an exponential solution: 

M stars (t) = M sc exp [(* - t c )M, (5) 

where the characteristic time r = 5~ l . Equation (5) remains valid during the supernovae 
phase of the parental cluster. Within this time interval, the continuous re-processing of the 
ejected material drives, as shown in Figure 4, the star formation rate to grow exponentially: 

SFR{t) = = ^ exp [(t - t c )/r\. (6) 

at r 
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Note that the characteristic time r is defined by the coeval star cluster model and is 
applicable to both the initial cluster and also to every new generation of stars formed from 
the re-ejected material. If one assumes that the energy and mass deposition rates remain 
constant during a star cluster life time and that v^, the resultant terminal speed in the 
absence of radiative cooling, equals 1500 km s _1 for all clusters, upon normalizing to the 
energy deposition rate adopted for the parental 10 7 M cluster (d = lO 41 Mj/lO 7 M erg 
s _1 ), then r = Mj/Mj = M i vl /(2C i ) = 7.5 x 10 7 yr, where Mj is the mass in each stellar 
generation, including the initial or parent cluster. The total mass in stars born within the 
star cluster volume, over the characteristic time r ~ 40 Myr, is given by equation (5). 



4. Conclusions 

We have here shown that massive and compact coeval clusters instead of driving a 
superwind able to disperse the surrounding ISM and even channel its way into the IGM, 
events that have make them been regarded as negative feedback agents, they are in fact 
extreme examples of positive star formation feedback. 

Massive and compact coeval clusters appear in the Csc vs size diagram above the 
threshold line, in the region where radiative cooling inhibits the development of stationary 
superwinds. We have here shown that in such cases the matter reinserted, through stellar 
winds and supernovae, is unable to escape and that after a short phase of matter accumu- 
lation, a new stationary solution in which M becomes equal to the SFR is rapidly met. A 
positive feedback condition in which new stellar generations result in situ, from the collapse 
of the matter reinserted by the sources evolving within the star cluster volume. 

The secondary star formation process while causing a faster mass deposition rate, drives 
the SFR to grow from 0.1 to 0.25 M Q yr -1 over the parent cluster supernova phase (~ 40 
Myr). The continuous re-processing of the ejected material leads effectively to a continuous 
transformation of the high mass stars into a low mass (< 8 M Q ) population, keeping con- 
stant the total mass of the stellar component. In this way massive clusters may survive as 
gravitationally bound systems (see discussion in Dale et al. 2005 and Martin- Hernandez et 
al. 2005 and references therein). 

A central issue, regarding ISM studies, is the fact that the more massive and compact 
clusters (as those detected by HST and large ground-based telescopes), are unable to generate 
superwinds and shed their metals into the ISM or the IGM. Their evolution leads to many 
stellar generations and thus to a mixture of stellar populations, all contaminated by the 
products from former stellar generations. Exacerbated episodes of star formation that leave 
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no trace of their evolution in the ISM. Clearly, if the most massive star formation events do 
not contribute to the chemical evolution of their parent galaxy, then several issues, such as 
the mass-metallicity relationship between galaxies, ought to be revised. 

We thank our anonymous referee for comments and suggestions. This study has been 
supported by CONACYT - Mexico, research grant 47534-F and AYA2004-08260-CO3-O1 
from the Spanish Consejo Superior de Investigaciones Cientificas. 
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Fig. 1. — The threshold line. The location of the threshold line, the line that divides the Csc 
vs size (Rsc) diagram into two distinct areas. There the matter deposited within the SSC 
volume, through winds and SNe, ends up either streaming away as a stationary (adiabatic or 
radiative) wind (below the line) or it accumulates to end up steadily being driven into new 
episodes of star formation (above the line). Two different locations of the threshold line, for 
an assumed full thermalization efficiency (e = 1), are display for an ejecta metallicity value 
equal to solar and ten times solar. The change in location for different values of e is also 
indicated. 

Fig. 2. — Positive star-forming feedback. The figure shows the rapid growth of the reinserted 
gas density within the SSC volume (rising solid lines), interrupted as it reaches the gravita- 
tional instability criterium (pj), leading to a new stellar generation and to a new phase of 
matter accumulation. The situation changes slightly as the available photon flux becomes 
unable to sustain all the SSC volume fully ionized, allowing for recombination and further 
cooling of the recombined matter, and thus to a much lower limit of the Jeans instability 
criterium. This allows the accumulation time to become comparable to the free-fall time, 
establishing a new stationary condition in which M = SFR. The density value {p gas ) required 
for this condition is also indicated in the figure as well as the maximum density value (phii) 
that can be supported fully ionized within the SSC volume, with the available UV photon 
flux. 

Fig. 3. — Figure (a) show the evolution of the UV photon flux produced by the main star 
cluster (dotted line) and the total flux resultant from the contributions from all further 
stellar generations (solid line), b) displays fnn, the fraction of the SSC volume that can be 
kept fully ionized with the available UV stellar photon output. 

Fig. 4. — Positive star-forming feedback. The evolution of the star formation rate throughout 
the life-time of the main cluster. 
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